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Abstract 


In the present thesis, we study the heat flow in mesoscopic one-dimensional transport sys¬ 
tems. Using the analysis of full counting statistics, we calculate the cumulant generating 
function of the particle and heat flows and prove its symmetry. The symmetry produces the 
relations among transport coefficients of the particle and heat flows when we expand these 
flows with respect to the appropriate affinities. Moreover, we consider the generalized flows 
which are superpositions of the particle and energy flows. We show that we can choose the 
appropriate affinities of the generalized flows and derive the relations among their transport 
coefficients when we expand the generalized flows with respect to their affinities. 
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Chapter 1 

Intoroduction and Notation 


In this chapter, we explain a brief history of studies on linear irreversible thermodynamics and 
its application to thermoelectric devices. We then introduce mesoscopic transport systems, in 
particular mesoscopic thermoelectric systems, which can be beyond linear-response regime. We 
hnally explain notations which are used throughout this thesis. 


1.1 Introduction 

1.1.1 Thermoelectric device as a heat engine and linear irreversible 
thermodynamics 

Thermoelectric devices, which convert heat to work or vice versa, have helped the development 
of our society. Applications of thermoelectric devices include thermoelectric generator, thermo¬ 
electric refrigerator, and so on [1] . In order to analyze thermoelectric devices, linear irreversible 
thermodynamics has been used HH. Linear irreversible thermodynamics is a phenomenological 
formalism which was constructed mainly by Onsager. Although the standard thermodynamics 
can treat only equilibrium, not transport problems, the linear irreversible thermodynamics lets 
us handle the latter. 

Let us overview Onsager’s formalism of linear irreversible thermodynamics (Ml. We con¬ 
sider two regions whose temperatures and chemical potentials we can dehne as shown in Fig. 11.11 
We then assume that the system has reached a non-equilibrium steady state in which there are 
constant flows from one region to the other. Let us expand the particle and energy flows, Jn 


temperature 

temperature 

T + AT Jn 

^ T 

chemical 

chemical 

potential h 

^ potential 

lJi + Afx 



Figure 1.1: The system in consideration. 
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and Je, in terms of yu/T and —1/T: 


Jn = LtvatA^^ j + Lne^ 
Je = LEN^iy^^ + Lee^ 



( 1 . 1 ) 


( 1 . 2 ) 


The expansion coefficients L^e, Len, and Lee are the transport coefficients for the 

particle and energy flows. The Onsager-Casimir relations are the relations among them under 
a magnetic held B [3HS]: 


Lnn{B) = Len{—B), 

Lne{B) = Len{—B), 
Lee{B) = Lee{—B). 


(1.3) 

(1.4) 

(1.5) 


Recently, the thermoelectric device has been studied as a heat engine from a point of view 
of the efficiency at the maximum power [SHig. Let us describe the setup to consider the 
thermoelectric device as a heat engine. We set the chemical potential of the right reservoir 
to be higher than the left, while the temperature of the left reservoir to be higher than the 
right so that an electric current may go from left to right against the difference of the chemical 
potential; in other words, we set AT > 0 and A/i < 0 in Fig. 11.11 What happens per unit time 
is the following. Electrons gain heat dehned by Jq = Je — lJn from the hot left reservoir, go 
to the right against the potential difference — A/r = |A/i|, during which electrons do the work 
of amount J]sf\Afi\. We can thus consider this system as a heat engine. The efficiency rj of this 


engine is therefore given by 


Jn\A^\ 


( 1 . 6 ) 


In order to analyze the efficiency of a heat engine, linear irreversible thermodynamics is still 


useful. When the relations for the particle and energy flows, Eqs. fll.3p - fll.5p . are valid, we can 
prove that the Onsager-Casimir relations are also valid for the coefficients of the particle and 
heat flows [2]. In order to prove it, we expand Jjy and Jq as follows: 



(1.7) 


and 


Jq — Je — lJn 












Let us define new transport coefficients as follows: 

Jat = + GatqA^-—^ , (1.9) 

+ GqqA^- —^ . (1-10) 

Comparing Eqs. fll.7p -- fll.8p and Eqs. fll.9p - fll.10p . we can express the new transport coefficients 
Gnni Gnq, Gqn, and Gqq in terms of Ljsiq, Lq^, and Lqq: 

Gnn = Lnn, ( 1 - 11 ) 

Gnq = Lme — ( 1 - 12 ) 

Gqn = Len — (1-13) 

Gqq = Lee — I^{Len + Lne) + I^^Leie!. (1-14) 

We can then prove the Onsager-Casimir relations for the particle and heat flows under a mag¬ 
netic held B as follows: 

Gnn{B) = Lein{B) = Le!n{—B) = Gein{.—B), (1-15) 

Gnq{B) = Le[e{B) — ^Lein{B) = Len{—B) — ^Le[n{—B) = Gqei{—B), (1-16) 

Gqq{B) = Lee{B) — ^{Len{B) + Leie{B)) + ^^Lein{B) 

= Lee{—B) — ^{Leie{—B) -|- Len{—B)) + ^^L]^ei{—B) = Gqq{—B), (1-17) 

where we used the Onsager-Casimir relations for the particle and energy hows, Eqs. fll.3p - fll.5p . 
We thus conclude that if the Onsager-Casimir relations for the particle and energy hows are 
valid, the Onsager-Casimir relations for the particle and heat hows are also valid. Using these 
relations, we can further prove that the upper limit of the efficiency Eq. fll.6l) is the Carnot 
efficiency, using the Onsager-Casimir relations and the positivity of the entropy production [7] . 

1.1.2 Thermoelectric device in mesoscopic transport systems 


Reservoir 


Reservoir 

Chemical Potential 


Chemical Potential 

P-L 

Temperature 

a quantum wire 

Ur 

Temperature 

7r 



Figure 1.2: The mesoscopic system which has a quasi one-dimensional wire and two reservoirs 
attached to it on both sides. 

Mesoscopic transport systems are systems in which a conductor of length Lg much shorter 
than the momentum-relaxation length Lp and the phase-relaxation length L^, is attached to 
reservoirs [TTlfT^ . Various interesting phenomena occur because of the length scale, which is 
called the ballistic transport regime. Let us consider here the simplest system shown in Fig. 11.21 
which consists of a quasi-one-dimensional quantum wire in the ballistic transport regime and 
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two reservoirs attached to it on both sides. The Landauer-Biittiker formula is particularly 
useful in this regime [I3MT6] : 



( 1 . 18 ) 


where I is the electric current across the system, Eq the ground-state energy of the wire, r(e) 
the total transmission probability of the wire, and /a(e) the Fermi distribution function of a 
Fermi gas in a reservoir a = L, R; see the notation section 1.2 and Chapter 2 for details. Such 
a system was theoretically considered by R. Landauer in 1957 [13], but had not been realized 
experimentally until 1988 d- Thanks to the improvement of sub-micron technology today, 
more rehned mesoscopic transport systems are made experimentally, which is also stimulating 
vigorous theoretical researches. 

With the development of the research in the mesoscopic transport systems, the mesoscopic 
thermoelectric device has also been considered theoretically and experimentally [91- [T^IT8] . It is 
expected to have a good efficiency because of its little heat leakage [T8] . 

In such a device, nonlinear effects can occur easily. For example, an experiment [19] suggests 
that the Onsager-Casimir relations, which are valid in the linear-response regime, are broken 
under a strong external held. This implies that one can make mesoscopic thermoelectric devices 
in a nonlinear regime under controlled external helds such as the difference of the chemical 
potential and the temperature. However, most of the theoretical approaches are still in the 
linear-response regime [Mg. There are a limited number of researches in the nonlinear regime 
[20ll2T| . but a general nonlinear theory for mesoscopic thermoelectric systems, which would be 
a counterpart of Onsager’s formalism in the linear-response regime, is yet to come. 

Indeed, Saito and Utsumi [22] have recently found relations among nonlinear transport 
coefhcients of the particle and energy hows, using full counting statistics [231121]. In the linear- 
response regime, if the Onsager-Casimir relations for the particle and energy hows are valid, 
those for the particle and heat hows are also valid as we showed above. There should thus be 
similar relations among nonlinear transport coefhcients of the particle and heat hows. We have 
indeed found them using full counting statistics, which we explain in this thesis. 

In Chapter 2, we review the Landauer-Biittiker formula [T3IIT6] . which is essential in treating 
mesoscopic transport systems. Using it, we calculate the average particle how as well as its 
second-order cumulant. In Chapter 3, we review the work by K. Saito and Y. Utsumi [22], in 
which the authors obtained the relations among the transport coefhcients of the particle and 
energy hows. Using the analysis of full counting statistics [231121], we calculate the cumulant 
generating function of the particle and energy hows and prove its symmetry. The symmetry 
produces the relations among transport coefhcients of the higher-order cumulants of the particle 
and energy hows. In Chapter 4, we properly dehne the heat how in mesoscopic transport systems 
and derive the relations among the transport coefhcients of the higher-order cumulants of the 
particle and heat hows. Moreover, we introduce the generalized hows which are superpositions 
of the particle and energy hows. We show that we can choose the appropriate affinities of the 
generalized hows and derive the relations among their transport coefficients when we expand 
the generalized hows with respect to their affinities. We hnally discuss their application to the 
calculation of the nonlinear Seebeck coefficient. In Chapter 5, we summarize our results and 
discuss possible future works. 
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1.2 Notation 


Let us here fix the notation. Throughout the present thesis, k-Q denotes the Boltzmann constant, 
h the Planck constant and e the elementary charge; for simplicity, we put = 1 throughout 
this thesis. 

Let fa{^) denote the Fermi distribution function of a Fermi gas in a reservoir a: 

= + (1.19) 

where (3a is the inverse temperature 1/T„ with Ta the temperature and /io, is the chemical 
potential. 
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Chapter 2 

Landauer-Biittiker Formalism 


In this chapter, we explain the Landauer-Biittiker formalism, calculating the first- and second- 
order cumulants, namely, the average and the noise of the particle flow. In Section 2.1, We 
present an elementary introduction of the Landauer-Biittiker formalism because we would like 
readers to understand the formalism intuitively hrst. In Section 2.2, we introduce the approach 
using the second quantization to calculate the average as well as the noise of the particle flow. 
Note that we neglect the spin of electrons throughout this chapter. 


2.1 Landauer-Biittiker formalism for calculation of av¬ 
erage current 


The Landauer-Biittiker formula is used to calculate the electric current through a quasi-one- 
dimensional mesoscopic conductor. The word ‘quasi-one-dimensionah indicates a conductor 
with the x-direction free while y- and 2 ;-directions conhned. We specihcally consider the system 
shown in Fig. 12.11 which consists of a quasi-one-dimensional quantum wire and two reservoirs 
attached to it on both sides. 

We here consider free electrons; that is, we neglect the electron-electron and electron-phonon 
interactions. The Schrodinger equation of an electron in the wire is given by 


H^{x,y,z) 


2m 


+ V{y,z) 


'^{x,y,z) = E^{x,y,z), 


( 2 . 1 ) 


where H is the Planck constant, k is the wave-number vector, m the effective mass of electrons, 
V {y, z) the conhned potential, and E the energy. Let us here separate the variables of \k(a:, y, z) 


Reservoir 


Reservoir 

Chemical Potential 


Chemical Potential 

P-L 

Temperature 

a quantum wire 

Temperature 

7r 



Figure 2.1: The system in consideration. 
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Figure 2.2: Dispersion relation in the x direction. 


as follows: 

=^p{x)cl){y,z), (2.2) 

which gives the Schrodinger equation in each direction in the forms 

(2.3) 

-^ ^ z) = Ey^,(j){y, z), (2.4) 

where each of and Ey^^ is the energy in the corresponding direction. Solving Eqs. fl2.3p and 
( El , we obtain the wave function '^{x,y,z) as follows: 

y, z) = z), (2.5) 

where is the wave number in the x direction given by 

^ ^2.6) 

Note that as the electrons are conhned in the y and 2 ; directions, they have the discrete energy 
Ey^z = -E'n, where n is the label of the level. Let us refer to the levels labeled by n as ‘channels’, 
through which electrons are transported in the x direction. The total energy of an electron is 

E = E^ + Ey^, = ^ (2.7) 

see Fig. 12.21 Throughout this thesis, we denote the energy for /c^. = 0 and n = 0 by the 
ground-state energy Eq. 

In the Landauer-Biittiker formalism, the important quantity is the transmission coefficient 
of the wire. We therefore explain the quantum scattering problem for a while, particularly the 
S-matrix. 
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Figure 2.3: The definition of the amplitudes a“, a°“*, 6'“, and 


2.1.1 Single-channel case 

Let us first consider the S-matrix in the case of a single channel, in which we denote a™, a°“*, 
and as the amplitudes of the incoming flow from the left reservoir, the outgoing flow to 
the left reservoir, the incoming flow from the right reservoir, and the outgoing flow to the right 
reservoir, respectively (Fig. 12.3p . The relation among the amplitudes a™, a°“*, 6™, and is 
expressed by the S-matrix in the form 



( 2 . 8 ) 


where r and r' are the reflection coefficients while t and t' are the transmission coefficients. We 
here remark that and |rp are the transmission and refiection probabilities from left to right 
and |Fp and |r'p are those from right to left. 

We can prove that the S-matrix is unitary as follows. The conservation of the flux gives the 
conditions 


I 2 


+ = 


.^out 1 2 


+ \b' 


out 1 2 


(2.9) 


Using Eqs. fl2.8p and fl2.9p . we have 


+ |6'“P = (a'“* t'"') = (a“‘t , (2.10) 


which gives = 1 and S~^ = = 1. The second condition gives = 1. The 

S-matrix is therefore unitary. 

Using these unitarity conditions, we can derive the following relation among the transmission 
and refiection coefficients: 


/r t'\ /rf _ f\r\^ + \t'\‘^ + _ fl 0\ 

\t r'j \t'^ ry v’ 


( 2 . 11 ) 


which gives the condition 


r|2 + |F|2 = |rf + |f|2 


The relation = 1 also gives the condition 


1 . 


rp + |tp = |rf+ |ff = 1. 


( 2 . 12 ) 

(2.13) 


Using these conditions, we obtain the relations 


r|2 = |rf and \t\^ = \tf 


(2.14) 
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Figure 2.4: The definition of the elements of a'“, and 6°“h 


2.1.2 Multi-channel case 

We can easily extend the single-channel case to the multi-channel case. Let N denote the 
number of channels. Note that N can be infinite in principle, but in real materials, the electron 
has a cut-off energy because of the energy band structure, and hence N may be finite. We at 
first define the iV-dimensional vectors a'°, h™ and whose components a'^, 
and 6°'^* for 1 < n < iV, respectively, denote the amplitudes of the incoming flow from the left 
reservoir, the outgoing flow to the left reservoir, the incoming flow from the right reservoir, and 
the outgoing flow to the right reservoir, each in the nth channel (Fig. 12.41) . We can express the 
relation among these vectors with the S-matrix in the form 



(2.15) 


where r and r' are N x iV-dimensional reflection matrices and t and t' are N x iV-dimensional 
transmission matrices. 


Similarlv to Eqs. f|2.9p-f|2.10p. we obtain the unitaritv condition SS^ = S = 

1, which 

gives the conditions of r, r', t, and f: 



N 

N 


SS' = l«^(|r.,|^ + |('yP 

) = = 1; 

(2.16) 

i=i 

i=i 


N 

N 



= Ed’-t.P + IPit) = 1- 

(2,17) 

i=i 

i=i 


They are followed by 



N 

N 


n = 


(2.18) 

i=i 

i=i 



where r^, r'ij, ttj, t'ij are components of the reflection and the transmission matrices r, r', t, 
t', respectively. We here remark that |r.yp and \r'ij\‘^ are the reflection probabilities from the 
jth channel of the left and right to the ith channel of the left and right, respectively, while 
and |FpP are the transmission probabilities from the jth channel of the left and right to 
the ith channel of the right and left, respectively. Summing both sides of Eq. fl2.18l) over i, we 
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obtain the equality 


N N 

i,j=l i,j=l 

N 

= 

*J=i 


(2.19) 


where in the second equality we replaced Vji and tji with and tij, respectively. We can obtain 
from Eq. fl2.19p the following equality of the transmission probabilities: 


N N 

Y, I4f = Z i*«it 

i,j=l *,J=1 


( 2 . 20 ) 


where the left-hand side is the transmission probability from left to right and the right-hand 
side is that from right to left. This equality is a generalization of Eq. fl2.14l) to the multi-channel 
case. Note that the transmission probability 02.201) can be expressed in the form 


N 


= Tr(ttt). 

*j=i 


( 2 . 21 ) 


2.1.3 Calculate the current 

Let us now calculate the electric current within the Landauer-Biittiker formula. We calculate 
it in the multi-channel case below. We make the following assumptions in order to do so: 

• The current coming into the lead holds the Fermi distribution of the reservoir in which 
it originally was and relaxes in the reservoir which it goes into. 

• The current which goes from the lead into the reservoir is not reflected back into the lead. 

• The electrons in the lead are one-dimensional non-interacting Fermi particles, and there¬ 
fore the current to the left and one to the right are independent of each other. 

Under these assumptions, the current which flows in the zth channel in the energy range [e, e+de] 
is given by 

d/f (e) = ev{e)Ue)D,{e)Tt{e)de, (2.22) 

where e is the elementary charge, a = L, R denotes the current from left to right and that from 
right to left, respectively, v{e) = de/d/c is the group velocity of electrons, Di{e) = dk/de = 
{hv{e))~^ is the density of states of one-dimensional ideal Fermi gas, and is the transmis¬ 

sion probability for electrons to transmit from the Rh channel in the left lead to a channel in 
the right lead, while the opposite. 

We can express T^{e) in terms of the elements tij of the transmission matrix; for example, 
we have 

N 

= p-23) 
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because the matrix element tji is the probability amplitude with which the incoming wave in 
the ith channel in the left lead transmits into the jth channel in the right lead as we explained 
below Eq. fl2.18p . We can similarly express Th(e) in terms of t'ji in the form 

N 

(2,24) 

i=i 

The total current in the energy range [e, e + de] is the difference between the left-going 
current and the right-going current: 

N N 

i=l i=l 

N 

= jE(A(f)fe(f)t-/R(f)l‘t(')l')* 

N 

= ^Efe(')l'(A(f)-/R(f))*. (2.25) 

where in the last equality we exchanged the dummy variables i and j and used Eq. fl2.20l) . We 
obtain the total current by integrating dl with respect to the energy e: 


1 = j dl{e) 

■ N 

El%(')t(A(')-A(')) 

.*J=1 

= t/ c^er(e)(/L(e) -/R(e)), (2.26) 

where r(e) = Yl!ij=i = Tr(ttl) is the total transmission probability. Equation fl2.26p is 

called the Landauer-Biittiker formula [T3HT6]. 



2.1.4 Conductance quantization 

Let us observe the conductance quantization within the Landauer-Biittiker formula at zero 
temperature. The Fermi distributions in the left and right reservoirs at zero temperature are 
respectively given by 

/„(e) = 0(e-p„), (2.27) 

where 0(e — fia) is the step function. The current I is given by the Landauer-Biittiker formula 
(I2.26p . which in the present case reduces to 

e 

I = - der(e)[ 0 (e - pl) - 0 (e -/^r)]- (2.28) 

" Jeg 

Let us calculate the conductance G, which is dehned by 
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(2.29) 
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where V is the voltage difference dehned by Id = (/il — fiR)/e. Using Eq. fl2.28p . we have 


G=^ f der(e) ^[0(e-/iL)-0(e-/iR)] 

= (2.30) 

where Ep = fip = f^R denotes the chemical potential of the left and right reservoirs at U = 0, 
namely the Fermi energy. 

Let us consider the transmission coefficient t{Ep). Because of the zero temperature, there is 
no electron which has an energy above Ep, so that no electrons transmit from and to a channel 
whose bottom energy is higher than Ep. This fact gives the transmission function 

n 

t(Ef) = Y. \t‘iiEF)?. (2.31) 

ti=i 

where n is the number of channels below Ep. Assuming that the diagonal elements \tii{Ep)\'^ are 
all equal to a constant T and the other elements are equal to zero, we observe the conductance 
quantization: 

G = —Tn. (2.32) 

h 

We can End from Eq. fl2.32p that the conductance increases stepwise by a discrete value of 
{e^/h)T as the number of channels below the Fermi energy increases. This conductance quan¬ 
tization was indeed observed in the experiment im. 

2.2 Calculation of the noise with second qnantization 

In this section, we calculate the shot noise within the Landauer-Biittiker formalism. In order to 
do this, we use the second-quantization approach [T6l[25]. First, we derive the average current, 
whose result is the same as the Landauer-Biittiker formula fl2.26p . We then derive the expression 
of the shot noise. The advantage of this approach is that the Pauli exclusion principle comes 
in naturally. 

2.2.1 Average current 


X = 

= 0 





< - bk 


lead 

scatterer 


Figure 2.5: The lead connected to the scatterer. We denote and a\ the annihilation and 
creation operators of the right-moving electrons, while bk and the annihilation and creation 
operators of the left-moving electrons. 
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Let us first find the current operator of the system shown in Fig. 12.51 The Hamiltonian of 
the lead is 


H = y^efcc|.Cfc, 
k 


(2.33) 


where Ck and c\ are the annihilation and creation operators, respectively. The energy of the 
electron is given by 




2m ’ 


(2.34) 


where k is the wave number of the electron and m the effective mass. We here remark that 
we may derive the Landauer-Biittiker formula without assuming the dispersion relation fl2.34p . 
but for simplicity, we assume that the dispersion relation is given by Eq. fl2.34p . The current 
operator is given by 


I{x) 


eh 

2mi 



dV'^ 

dx 


'ip{x) 


eh 

2imL 


k^k' 


eh 

2mL 


k,k' 


where ijj{x) and 'ijj^{x) are the held operators given by 


^ k 


(2.35) 


(2.36) 

(2.37) 


Since the current is conserved, Eq. fl2.35p should be independent of x, and hence should reduce 
to 


Hx) = m 


eh 

2mL 


k,k' 


(2.38) 


In the Landauer-Biittiker formalism, we consider the current going to right and one going 
to left independently, which enables us to dehne the annihilation and creation operators of the 
right-moving electrons, and a|,, as well as the annihilation and creation operators of the 
left-moving electrons, and 6^. Using these operators, we can express the Hamiltonian in the 
form 


H = Y ^kaittk + Y ^kblbk- 

k>0 k<0 


(2.39) 
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a: = 0 

I 


OLfc-► 

bLk< - 

leadL 


scatterer 


*-Qr/c 

► bRk 
leadR 


Figure 2.6: Two leads attached to the both sides of the scatterer. 


We can write down the current fl2.38p in terms of the new operators as 


eh 


I{x) = /(O) = + k')clck' 


k,k' 


eh 

2mL 


{k + k')clck' + {k + k’)c\ck' + {k + k’)c\ck' 

.k,k'>0 k^k'<Q k^k'&X.kk' <0 


e 

2L 

e 

TL 

e 

TL 

e 

2L 


(wfc+yfc')4cfc' + (^fc+wfcOcicfc' + 

_k,k'>0 k,k'<0 k,k's.t.kk' <0 

y^ (wfc+ffc')(4cfc' - cyc_fc') + (ufc - ffcOcic-fc' + {vk' - Vk)c^_f.ck' 


k,k'>0 


y^ (Wfc + ~ blbk>) + {vk - Vk')albk' + {vk' - 


k,k'>0 


{vk + Vk'){alak> - blbk>) + (cross terms), 

k,k'>0 


(2.40) 


where Vk = hk/m is the group velocity of the electrons. Note that the cross terms vanish when 
k = k'. 

We next consider the system where two leads are attached to the both sides of the scatterer as 
shown in Fig. 12.61 We denote the annihilation and creation operators of the incoming electrons 
in the lead a = L,R by aa,k and respectively, while those of the outgoing electrons by ba,k 

and 4,fc- 

The S-matrix dehned in Eq. fl2.8p gives the relation 


( I'-A = s (“■'d, 

\bR,k J \aK,k J 

where we now use the notation 

c.^/^^LL(fc) SMk)\^fr f\ 
V^RL(fc) SnK{k)J \t r'J- 


(2.41) 


(2.42) 


The unitarity condition of the S-matrix is guaranteed by the following commutation relations: 


,0'a,k, 

^ol' ^k^ ^OL^OL'^k^k' ^ 

(2.43) 

P‘a,k^ Old',k'\ 

,k'\ 0, 

(2.44) 



(2.45) 
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We can express the current operator Ii^{x = 0) going from the lead L to the scatterer as 
Eq. (ICTD : 

/l(x = 0) = ^ ^ {Vk + 'yfcO(«L,fc®L,fc' - + (cross terms) 

k,k'>0 

= +'^^0 
k,k'>0 

t 


X 


®L,fc®L,fc' (*S'LL(fc) flpfc + 5 'lr(A;) )®L,A;' + *S'LR(fc )(^li^k' 

- (cross terms) 

('^^fc + nfc')«I,fc"4L^(^> /;;')«/?,fc' + (cross terms), 








2L 


(2.46) 


q=L,R /3=L,R k,k'>0 


where 


Af{k, k') = 4,a5L,/3 - S^kySi^pik'). (2.47) 

Let us introduce the statistical average of the creation and annihilation operators: 

{^a,k^0tk'') ^a,)3^k,k' fai,k'), (2.48) 

which means that the distribution of incoming electrons is the Fermi distribution of the bath 
where they originally were. We thereby calculate the statistical average of the current as 

{h) = y] y y ('^fc + Vk')Al^{k, /c')(«I,fca/3,fc') + (cross terms) 


a=L,R /3=L,R k,k'>0 


^ E Y.^vtAr(k,k)uk) 

q;=L,R k>0 

- de ^r(e,e)/a(e), 

" JEa „±rD 


(2.49) 


q=L,R 


where in the second equality the cross terms vanish because they survive only when k y k' but 
the statistical average insists k = k'. In the last equality we replaced the summation with 
the integral with respect to energy, assuming that the value of L is large enough for us to do 
the transformation as follows: 


where 


E 


D<(e) = 


-)■ 


dk 


1 dk 


deDi{e), 

(2.50) 

m 

u 2'Kh^k 

(2.51) 


is the density of states of the one-dimensional ideal Fermi gas. 

Using the relations between the elements of the S-matrix, we have 

/lp(£,€) = 1 - Sll(£)-Sll(£) = 1 - |r|^ = = t(£), 

yl!“(e,6) = -Si,R(e)*SLE(6) = -|if = -|«p = -r(e), 

which reduces (12.491) to the Landauer-Biittiker formula 


(2.52) 

(2.53) 


(4) = 


h 


der(e)(/L(e) - fnie)). 


Eg 


(2.64) 
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2.2.2 Shot noise 


The current noise contains much important information of the system in question. The power 
spectrum of the current noise is dehned by 

T T 

5(a;) = lim |;|A/(a;)p = lim / dt [ (2.55) 

T^oo I T^oo I T T 

''-2 -^-2 

where u is the frequency, T is the time to observe the current, I(t) is the electric current, 
A/(t) = (/(t)) — I{t) is the fluctuation of the electric current from its average value, and 
I{oj) = f%dtl(t)e‘^ is the Fourier transform of the electric current To calculate the 

shot noise, let us dehne the correlation function of the electric current as follows: 


C{t,t') = {Ah{t)Ah{t')) = {h{t)lL{t')) - (hit)){hit')), 


(2.56) 


where A/L(t) = hit) — (hit)) is the fluctuation operator of the current. Using this function, 
we can express the power spectrum of the current noise, Eq. (I2.55p . as follows: 


Siu) = lim — 
^ ^ T^oo T 


T 
■ 2 


dt / dt'Cit,t')e'‘^^'-''\ 


T 

2 


(2,57) 


which is of the same form as the classical noise power. Note, however, that the current operators 
hit) and hh) of different time do not commute with each other in the quantum case. When 
the Hamiltonian of the system does not depend on time, the system has the time-translational 
symmetry, and hence we can express the correlation function as Cit, t') = Cit — t'). Moreover, 
we assume that the correlation function Cit — t') tends to zero when the time difference \t — t'\ 
goes to inhnity. Based on these considerations, we can express the noise power Siuj) in the 
form: 


Siu) = lim — 


= lim — 

T-)-oo T 


~ lim — 

T^oo T 


, T 
2 


dt 


_T 

2 

T 

: 
2 
T 
' 2 


dt' 


dt' 


i-T 

r+oo 


T 
■ 2 


dt'Cit - 
dAtCiAt)e'‘^^' 
dAtCiAt)e'^'^' 


' —OO 


= 2 


dtCit)e 


iujt 


(2.58) 


In general, Siu) contains many components of different u. We, however, only calculate the 
zero-frequency component of «S'(a;) for simplicity: 


F(0) = 2 / dt(A/L(t)A/L(0)) = 2 / dti{hit)hm - {hit)){him ■ 


(2.59) 


In order to calculate 5'(0), we need the time evolution of the current operator hit)- We thus 
use the time evolution of the creation and annihilation operators: 


aa,kit) = e 

aUit) = 


(2.60) 

(2.61) 
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which gives us the expression of the current operator at time t: 


a=L,R ^=L,R k^k'>0 


where we ignore the cross terms because they will vanish when we take the statistical average 
below. Substituting Eq. 02.461) and its time evolution lL(t) into Eq. 02.59p . we can express 5'(0) 
in the form: 




dt ^ ^ (vk + Vk'){vk" + Vk'")Al^{k, k'") 

k,k',k",k'">0a,l3,a>l3' 


X 






(2.63) 


In order to calculate the statistical average of the creation and annihilation operators, we 
use Wick’s theorem |26j : 


(ABCD) = {AB){CD) - {AC){BD) + {AD){BC), (2.64) 


where A, B, C, and D are arbitrary Fermion operators. Substituting for A, ap^k' for B, 
aj^, for C, and for D, we obtain 






(2.65) 


In the second equality, we used the equation 




( 2 . 66 ) 


because the operators and ap^k'^p',k"' do not conserve the number of particles, so that 

the statistical averages of these operators become zero. Using Eqs. fl2.48l) and fl2.65p . we obtain 
the following relation; 


^a,/3' ^k,k'" ^/3,a'dk' ,k" fai^k'ji^l- fj^^ky^. (2.67) 


Using Eq. fl2.67p and the transformation fl2.50l) . we finally obtain the expression of the shot 
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noise as follows: 


^ ^ (^) / ^ ^ + Vk'){vk" + Vk'")Al^{k, k')Al^ (fc", fc'") 

^i{ek-e'^.)tlh 


X 


k,k',k",k"'>Qa,l3,a'P' 

J „ J „ \ iJ ^ mJ 


2 /*oo 

n^) E E (Wfc + + Wfc''')^L^(^) k')A1 ^ (fc", k'") 

k,k',k",k"'>0a,l3,a'l3> 


= 2 


e 


2L 

2 /*oo 



f Eq 


fc,fc'>0 a,/3, 

/*oo poo 

deDi[e) / de'Di[e') / dt 
J Eq J —oo 


X (6,6')xlf (6',e)/„(e)(l - 

a.,13 

0^2 ^oo 

/ *5;/lf(£,£)/lJ“(£,£)/„(£)(l-/4£)), 


( 2 . 68 ) 


where we used the relation 


^^^i{e-e')tlh 


27rh6{e — e') 


(2.69) 


in the last equality. 

In order to express the noise more simply, we use Eqs. fl2.52p and (12.531) as well as the 
relations 


e) = -SLL(e)'SLR(e) = -r’i', (2.70) 

7lp(e, e) = -SLR(e)-SLL(e) = -i"r, (2.71) 

which give the relations 

Al\e,e)Al\e,e) = Al\e,e)Al\e,e) = \r\^\tf = \r\^\tf = r(6)(l - r(e)). (2.72) 

We then calculate the integrand of Eq. fl2.68p as 
^ 7lL^(e, e)7l^"(e, e)/„(e)(l -/^(e)) 

a,/3 

= Al\e,e)A]:\e,e)h{e){l - h{e)) + e)7l^L(e, e)/L(6)(l - Ue)) 

+ /I^^(6,e)x4[:^(6,e)/R(6)(l - /L(e)) + 7l^^(6, 6)x4r(6, 6)/R(e)(l - fM) 

= ^(e)^[/L(e)(l - k{e)) + /R(e)(l - /R(e))] 

+ ^(e)(l - r(e))[/L(e)(l - f^e)) + /R(e)(l - /L(e))] 

= r(6)(/L(e) + Ue) - 2k{e)Me)) - r{enk{e) - /R(e))^ (2.73) 

which gives the hnal result of the expression of the shot noise: 

4p2 roo 

^(^) = -r [^(e)(/L(e) + fnie) - 2k{e)fii{e)) - r{ef{k{e) - fnie))^] . (2.74) 

IT' JEq 
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Chapter 3 

Full Counting Statistics 


In the previous chapter, we considered the hrst- and second-order cumulants of the particle flow. 
In this chapter, we consider the higher-order cumulants of the particle and energy flows. We 
review derivation [22] of the generalized Onsager relations among their transport coefficients 
using full counting statistics [23l|2l|. In Section 3.1, we introduce the full counting statistics 
[2^12^ briefly. In Section 3.2, we review Ref. [22]. We calculate the cumulant generating 
function of the particle and energy flows and prove its symmetry. This symmetry gives the 
generalized Onsager relations for the particle and energy flows. 


3.1 What is counting statistics? 


Counting statistics is an analysis method in which one counts a physical quantity in time s and 
examines its statistics. We would like to apply counting statistics to the transport phenomena 
driven by thermodynamic power (affinity), such as electrical conduction and heat conduction. 
We here explain a concept of full counting statistics [23l|2l] using the example of electrical 
conduction. 

Consider two reservoirs which have different chemical potentials /il and pr. The difference 
of chemical potentials, pl — /^r, causes an electric current. We here assume pl > /^r so that 
electrons may flow from right to left. We measure the number of the electrons Q for the time 
s. After we repeat the measurement many times, we obtain the probability distribution P{Q) 
of the number of electrons Q. Our aim is to know the probability distribution P{Q) for large 
enough s because we would like to know the non-equilibrium steady state, a long-time behavior 
of the system. 

We then explain how we obtain information of the cumulant from P{Q). Let us define the 
characteristic function Z{x), the Fourier transform of P{Q)'- 


Z{x)s^P(Q)e‘^'^, 

Q 


(3.1) 


where x is a variable called ‘counting field’. We can compute the expectation value of from 
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Z(x) as follows: 


d-Z(x) 


d(ixy 


d(ixy 




,*xQ 


x=o Q 

Q 

= m- 

Let us define the cumulant ((Q”)) of Q as follows: 

dr\ogZ{x) 


X=0 


((Q")) = 


diixy 


(3.2) 


(3.3) 


x=o 


For example, ((Q)) = {Q) is the expectation value of Q and ((Q^)) = (Q^) — (Q)^ is the variance 
of Q. 

Let us also dehne the cumulant generating function for the electric current: 


Fix) = ~ loS 

s—>-oo S 

We can, for example, calculate the cumulant of the electric current I as fallows: 


((/)) = 


dFix) 


diix) 


x=o 


,. (Q) 

= hm -, 

s—>-oo S 


which is the expectation value of the electric current, and 


an = 


d^Fix) 


diixy 


x=o 


= hm 

S^OO 


- (QY 


(3.4) 


(3.5) 


(3.6) 


which is the noise of the electric current. We thus obtain higher-order cumulants from Fix). 
Using this function Fix), we can obtain an asymptotic expression of PiQ) as follows: 


P{Q) = 


—oo 
roQ 


2tt 

1 


= hm — 

s—>-00 2 tI 


dxZix)e-^^^ 

dxy°sZ{x)-ixQ 

poo 


(3.7) 


where h(x) = Fix) — ixd and q = Q/s. We can evaluate the integral by the saddle-point 
method. Let x* denote the value of x which makes h(x) maximum: 


dhix) 


dx 


= 0, or 


dFix) 


x=x 


dx 


x=x" 


s 


(3.8) 


Expanding h(x) around X = X* to the second order, we arrive at an asymptotic expression of 
PiQ) as follows: 


PiQ) = hm ^ 

s—>-oo ZTT 


dx exp 

XD 

= lim = lim 


* I MV) + )x - x'Y ^ 


X=X 


(3.9) 






















3.2 


Review of ref. 



3.2.1 Setup and model Hamiltonian 



Figure 3.1: The system in consideration. 


We consider the model in which n quantum dots are connected to m reservoirs as shown in 
Fig. 13.11 The Hamiltonian of the model is given by 

m 

H = Hr+ H, + + Ht, (3.10) 

r=l 

where Hr denotes the Hamiltonian of the rth reservoir, the Hamiltonian of the quantum 
dots, Hint the Coulomb interaction between electrons on each dot as well as between the dots, 
and Ht the tunneling Hamiltonian: 

Hr = E (3.11) 

k(j 

where k is the wave number, a denotes spin, erk is the energy spectrum of the rth reservoir, 
®rfco- creation operator of the reservoir, and arka the annihilation operator of the reservoir; 


Hci = 


ia ija 


(3.12) 


where is the energy level of the ith dot, Uj is the hopping matrix between the dots, d\^ the 
creation operator of the ith dot, and dia the annihilation operator of the ith dot; 


Hint „ y ^ Uifxja'dj^^dj^idjijidifri 


IJCT(J' 


(3.13) 


where Ui^ja' is the strength of Coulomb interaction on each dot and between the dots; 


Hx ^ trkjdj^^Clrka H-C., 


rkicT 


(3.14) 
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where trki is the tunneling matrix between the rth reservoir and the ith dot. 

We later consider the situation where we measure the system at time t = —s/2 and t = s/2. 
We assume that the initial density matrix po at the time t = — s/2 is 


Po = Pd® Po 

_ exp[ 

Tl{exp[—l 3 r{Hr — PrNr)]}^ 


(3.15) 


where pd and p'^ are the density matrices of the quantum dots and the reservoirs, respectively, 
r the index denoting the reservoirs and d the quantum dots, (dr the inverse temperature of the 
rth reservoir, pr its chemical potential, and W its number operator: W = (^Ika^rka- We 
assume that there had been no interaction between the reservoirs and the quantum dots until 
the initial time t = —s/2 and hence their density operators are commutative then. 

We here assume that pd = l(i/2"', where Id is the identity matrix. Because we expect that 
the stationary state in the long-time limit is independent of the initial state of the dots, we can 
arbitrarily choose the initial state of the dots, and thus took pd = 1^/2"' for our convenience. The 
assumption becomes useful in proving the symmetry of the characteristic function Eq. fl3.3ip . 
whose details are shown in Appendix A. 

We then dehne the particle current operator and the energy current operator with the 
Heisenberg equation: 


dNr dd % [A^ f. , ^ ^ '^irkidj^^drka TH.C., (3.16) 

ik 

JH-p % \H-jp , -^t] ^ ^ '^^rk^rkidj^fjClrko' “hH.C.. (3.IT) 

ik 

Let us define the charge qnt and the energy qet of the rth reservoir during the measuring time 
s as follows: 


«/2 


QNr / dtJNr{t), 

( 3 . 18 ) 

J-sl2 


rs/2 


QEr = / dtJErit)- 

( 3 . 19 ) 

J -sl2 



3.2.2 Protocol of the measurement 

In order to obtain the expressions of the characteristic function Z{x) and the cumulant gen¬ 
erating function F{x), we need to know the probability distribution P{Q) as we discussed in 
Sec. 3.1. We can obtain the expression of P{Q), considering the following protocol of measure¬ 
ment [27] . 

1. At the initial time t = —s/2, we carry out a projection measurement on a reservoir r, and 
thereby obtain the result and the energy Q^riy of particles in it, which collapses 

the wave function of each reservoir into an eigenfunction of the particle and energy 

operators: 

Nr \^lC) = QnIu > 

Hr = QS. • 
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(3.20) 

(3.21) 





Note that the set of kets ^ complete set of the eigenspace of Nr and Hr. 

We make the same observation for the other reservoirs and obtain the set of particles 
{QNrv}T=i energies which collapses the wave function of the whole system 

into 

m 

|l'f‘> = 0IC'")®W“). (3.22) 

r=l 

where is the wave function of the quantum dots after measurement of all reservoirs. 


2. We let the whole system, which consists of m reservoirs and n quantum dots, evolve over 
time s with the unitary operator 

3. At the time t = s/2, we again observe a reservoir r to measure the number Q^r\ 

the energy Qj/l-x of the particles in it, which collapses the wave function of each reservoir 
into the eigenfunction I'l/r’x) of the particle and energy operators: 

Nr ICa) = Q%lx Ka) , (3.23) 

Hr = Q%.,x Ka) • (3.24) 

Note that the set of kets {iV'r a)}^o i® ^ complete set of the eigenspace of W and 
Hr- We make the same observation for the other reservoirs and obtain the set of particles 
{Q^r x}T=i energies {Q%. x\T=ii which collapses the wave function of the whole system 
into 

m 

I'I'f > = 0 ® IV-J”), (3.26) 

r=l 

where is fhe wave function of the quantum dots after measurement of all reservoirs. 
After the measurement, we obtain the conditional probability Pu^xiiQNr}, {QEr}) of 
hnding the changes in the number QNr = QNr,u—QNr,x and the energy QEr = QEr,u—QEr,x 
of the particles in each reservoir under the condition that the initial state of the whole 
system is hxed to 


Hu^x(^{Q Nr }) {QEr} ) 




Kd'f 


Z>',A S.t. QjSlr — Qnv^u QNr,X^QEl - Qet^u QEr,\ 

2 


A=0 


X - («“„ - Q«r,A)]3|QE. - (Qfl - Q|".a)]. (3.26) 

r=l 


4. We iterate the protocols 1 to 3 for the same initial density matrix po and obtain the 
probability of hnding the changes in the number Q^r and the energy Q^r of the particles 
in each reservoir with the initial state obeying the grand-canonical distribution 

(4/™* I Pol 4/™*) ill fhe form; 


^’(w».}.{Qi;j)= a^A({QA..},WBr})(>tr’iPoi>i'r*>. (3.27) 
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3.2.3 Symmetry of the cumulant generating function 

We can calculate the characteristic function Z(x) from the probability P{{QNr}, {Qet}) as 
follows: 


mXcrj, {XhrY B)= Y. P^iQNr). {Q Er}) J] (3.28) 

{QATrlifQBr} '”=1 

= , (3.29) 


where 


m 

V = JJexp 

r=l 


^i^Xcr^r P XhrHf) 
2 


(3.30) 


We show the details of the calculation in Appendix A. This characteristic function has the 
following symmetry 


^({Xcr}, {Xhr}] B) = Z{{-Xcr + iAcr}, {-Xhr + iAhr}', -B), (3.31) 

whose proof is also shown in Appendix A. 

Using this characteristic function fl3.29p . we can define the cumulant generating function of 
this model in the form: 

F{{Xcr},{Xhr}]B) = lim -\nZ{{xcT},{Xhr]]B). (3.32) 

s—>-cx> S 

The symmetry of the characteristic function (13.3Ih gives the symmetry of this cumulant gen¬ 
erating function as follows: 


^({Xcr}, {Xhr}] B) = F{{-Xcr + i^cr], {-Xhr + iAhr}] -B). (3.33) 

Noting that the cumulant generating function only depends on the difference between two 
counting fields [ 22 ] as in 


F{,{Xcr}, {Xhv}] B) = F{{Xcr}, B) , (3.34) 

where Xcr = Xcr — Xcr' and Xhr = Xhr — Xhr' with r’ hxed and r' r, we obtain the following 
symmetry of the cumulant generating function from the symmetry fl3.33p : 

F{{Xcr}, {Xhv}] B) = F{{-Xar + {-Xhr + iAhr}', -B), (3.35) 

where Acr and A^r are the affinities (thermodynamic forces) of the particle flow and the energy 
flow, respectively: 


Acr — |3r^^r — Pr'l^r'■> (3.36) 

Ahr = -^r + Pr', (3.37) 

with r' fixed to one of the affinities and r 7 ^ r'. We will explain why we choose these affinities 
in Chap. 4. The symmetry fl3.35p produces many interesting relations of transport coefficients. 
We next see the relations for the simplest two-terminal case. 
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3.2.4 Generalized Onsager relation in the case of two terminals 


Let us consider the case of two terminals, in which we can write the symmetry of the cumulant 
generating function in the form: 


F{Xc, Xh] B) = F{-Xc + iA, -Xh + iAh] -B), (3.38) 

where Xc = XcL - XcR, Xh = Xhh - XhR, A = A^l - A^r = - Af^R, and Ah = Ath - Ahr = 

—A + A- We can compute the kth cumulant of the particle flow Jjv and the energy flow Je 
by differentiating the cumulant generating function with respect to Xc and Xh as follows: 




F{xc,XhA) 

d{iXcY^d{iXhY^ 


X<^=Xh=Ac=Ah =0 


(3.39) 


We then dehne the transport coefficients as follows: 

,hM,n 


(3.40) 


For later use, we symmetrize and antisymmetrize the transport coefficients and the cumulant 
generating function with respect to the magnetic held: 


F±{Xc, Xh] B) = F(xc, Xh] B) ± F(xc, Xh] -B), 

rfcl,fc 2 / _ rhl,k2{JD\ _L rklMf ]D\ 

Fei,e2±\B) — ± ( F), 

which satisfy the relation 

d^^+>‘^+^^+^^FAXc,Xh]B) 


jkiM 


{B) = 


d{ixAA{%XhYAAcAA^ ^^=^,=^^=^,=0 
The symmetry in Eq. fl3.38p gives the symmetry of F±{xci Xh] B) in the form: 

F±{Xc, Xh] B) = ±F±(-Xc + iAc, -Xh + iAn] B). 


(3.41) 

(3.42) 


(3.43) 


(3.44) 


Partially differentiating both sides of Eq. fl3.44l) with respect to Xc times, Xh ^2 times, 
Ac l\ times, and Ah f '2 times, we arrive at the relation of the transport coefficients as follows: 


dkMiA+i2F^^^^^Xh]B) 

d{ixcfA{ixhAdAiAA^ 

_ ^ dk,+k2+i^+i2F^^_^^ + ^ ,Ah] B) 

d{iXcY^d{iXhY^dAi^dA^^ ’ 

which produces 

k\i:±(B) = ± E E Cn) (^')(-i)”‘*”=*‘‘'"‘"dTT';-t."±(B)' 

ni= 0 n 2=0 V 1 / V 2 / 


(3.45) 


(3.46) 


Note that the cumulant generating function F{xc, Xh] B) also depends on the affinities Ac and 
Ah when one differentiates F{xc,Xh] B). The relation fl3.46l) is among the coefficients with a 
hxed value of iV = fci + fc 2 + + £ 2 - For N = 2, we have 

L%{B) = Ll'ii-B), 

LlliB) = L\^,i-B), 

A,\{B) = Llli-B), 
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(3.47) 

(3.48) 

(3.49) 

















which are known as the Onsager-Casimir relations, as well as 


= 2L1S(B), 
iSS(S)=2L|>!(S). 


which are known as the fluctuation-dissipation theorem. For N = 3, we have 


L 

L 

L 

L 
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L 
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10 ,+ 
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01 ,+ 

30 

00 ,+ 


10 

20 ,- 


01 
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20 
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02 

10 ,- 


(B) 

(B) 

(B) 

(B) 

(B) 

(B) 

(B) 

(B) 

(B) 


rlO 

-^20,+ 

(B), 



rOl 

-^20,+ 

iB) = 

21“,+(B) - 

-LlUB); 

rlO 

-*^02,+ 

{B) = 

21”,+(S) - 

-lToAb), 

rOl 

-^02,+ 

{B), 



r03 

-^00,+ 

(B), 




Lll-jB) 

L^ol-jB) 

3 


6 






Lll_{B) + 2L'il_iB). 


(3.50) 

(3.51) 


(3.52) 

(3.53) 

(3.54) 

(3.55) 

(3.56) 

(3.57) 

(3.58) 

(3.59) 

(3.60) 


The relations among the transport coefficients of the cumulant of the particle flow 03.481) . O3.50p . 
03.52p . and 03.571) have been observed experimentally [281129] . 
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Chapter 4 

The Generalized Onsager Relations 
between Heat Flow and Particle Flow 


In the previous chapter, we consider the higher-order cumulants of the particle and energy flows. 
In this chapter, we treat the higher-order cumulants of the heat flow. We show the appropriate 
dehnition of the heat flow in mesoscopic transport systems and derive the generalized Onsager 
relations for the particle and heat flows. In Section 4.1, we review the several dehnitions of the 
heat flows, which might have confused researches in the past. We then explain the heat flow 
which seems to be the most relevant one. In Section 4.2, we hrst explain the affinities. We then 
derive the generalized Onsager relations for the particle and heat flows when we expand these 
flows with respect to the appropriate affinities. We also consider the generalized flows which are 
superpositions of the particle and energy flows. We show that we can choose the appropriate 
affinities of the generalized flows and derive the relations among transport coefficients of the 
flows when we expand them with respect to their affinities. For simplicity, we treat only the 
setup with two reservoirs throughout this chapter. 


4.1 Heat flow in mesoscopic transport systems 

4.1.1 Several ‘heat’ flows in mesoscopic transport systems 

We here introduce three expressions of the ‘heat’ flow in mesoscopic transport systems. The 
thermodynamic dehnition of heat is clear, but the wording of the ‘heat how’ may be confusing 
because there seems to be several dehnitions. 

The hrst dehnition of the heat how is given by 

Je = ]- [ der(e)e(/L(e) - /R(e)), (4.1) 

" JEg 

which is used widely in researches of mesoscopic transport [301132] . Although this ‘heat’ how 
should be called an ‘energy’ how, it is called a ‘heat’ how probably because it is considered in 
the situation where the energy how does not do work and hence all energy becomes heat. 

The second dehnition of the heat how is given by 

^ -/R(e)), (4.2) 

" JEg 
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which was often used in the dawn of the research of heat flow in mesoscopic systems [33] . This 
definition may have been taken from an equation in Callen’s textbook [2], 


Jq — Je — IJ'Jn) (4-3) 

where Je is the energy flow given by Eq. fl4.ll] . while J^r is the particle flow, for which the 
Landauer-Biittiker formalism gives the expression 

Jn = j- f der(e)(/L(e) - fnie)). (4.4) 

" JEg 

Substituting Eq. (14.11) and Eq. (14.4p into Eq. (14.31) . we would obtain the second definition (14.2p . 

Since this ‘heat’ flow was not microscopically derived, we do not know clearly where it 
flows. It is indeed ambiguous of which part of the system in Fig. 11.21 the chemical potential 
fi of Eq. fl4.3p is. We should probably choose fi so that Jq may satisfy Onsager’s reciprocal 
theorem. For example, in Ref. [33], the authors chose /i as (/il + Pr)/2 and in Ref. [33], the 
author chose /i as pl. The choices do not make difference in the linear response of the voltage 
difference (/tr — but differ in higher orders. 

The third definition is given by 

= T / deT{e){e - /ia)(/L(e) - fnie)), (4.5) 

E Jeg 

where « = L, R. This appears to be the same as the definition (14.4p with the chemical potential 
arbitrarily fixed, but we here make distinction because the definition fl4.5l] clearly specifies 
where the heat flow exists; we will show below that this definition gives the heat flow out of 
the left reservoir or into the right reservoir. 

4.1.2 Definition of the heat flow 

In this section, we derive the third definition of the heat flow (14.51) using the Landauer-Biittiker 
formula and thermodynamics. We here derive Jq specifically, but we can derive Jq in the same 
way. 


Reservoir 


Chemical Potential 
Ul 

quantum wire 

Temperature 

Tl 



Figure 4.1: The wire and the left reservoir. 

Consider the quantum wire with the left reservoir (Fig. I4.ip . We assume that the reservoir 
is so large that it is always in equilibrium. Hence we can define thermodynamic quantities of 
the reservoir such as the temperature and the chemical potential. 
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Since each reservoir is coupled only with the wire, all the heat generated in the left reservoir 
goes into the wire. Therefore, 

Ql = (4.6) 

where Ql is the heat generated in the left reservoir per unit time. The negative sign appears 
because we dehne the positive direction so that flows going to the right may be positive. 

Using the hrst law of thermodynamics, dQ = dU + dPU, we obtain 

Ql = Ul- Wl, (4.7) 

where Ul is the energy going in the left reservoir per unit time and Wl is the work per unit 
time done on the left reservoir. Let us use the framework in which the volume of the reservoir 
is constant. Hence, the pressure dose not do work on the reservoir. Therefore, 

Wl = MlJVl, (4.8) 

where Nl is the number of particles going in the left reservoir per unit time. 

In order to hnd Ul and Nl microscopically, we use the Landauer-Biittiker formalism, which 
gives 

^L = -j- [ der(e)e(/L(e) - fnie)), (4.9) 

n Jeg 

Nl = -j- f der(e)(/L(e) - fnie)). (4.10) 

n- Jeg 

We thus arrive at 

’^Q = T deT{e){e - /iL)(/L(e) - /R(e)). (4.11) 

" JEg 

We can derive Jq in the same way, expect that we have Qr, = Jq and therefore the sign of Ur 
and IVr are reversed, which results in 

1 

deT{e){e - /iR)(/L(e) - /R(e)). (4.12) 

n Jeg 

These heat flows Jq and Jq have not been observed experimentally as far as we know. In 
order to observe them, the experimental condition should be at low temperatures as was when 
the energy flow was observed [33]. We therefore show in appendix B the expansion of Jq and 
Jq at low temperatures. 

4.1.3 Constructing heat engine and its efficiency with the heat flow 

The heat flows Jq and Jq in Eq. 04.51) can be used to analyze mesoscopic heat engines [T2| in 
Fig. 14.21 We set the chemical potential of the right reservoir higher than the left, while the 
temperature of the left reservoir higher than the right so that an electric current may go from 
left to right against the difference of the chemical potential. What happens per unit time is 
the following. Electrons gain heat Jq from the hot left reservoir, go to the right against the 
potential difference, during which electrons do the work of amount IV, where I is the electric 
current and V is the voltage difference, and then dump heat Jq to the cold right reservoir. We 
can thus consider this system as a heat engine. Its efficiency rj is thereby given by 

IV 

V = 4 - (4-13) 
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Figure 4.2: Schematic picture of a mesoscopic heat engine. Note that we cannot dehne the chem¬ 
ical potential in the wire as this illustration might suggest because it is highly non-equilibrium 
there and therefore the thermodynamic quantities, such as the chemical potential cannot be 
dehned. 

The current I is given by the Landauer-Biittiker formula 

e 

I = -r c?er(e)(/L(e)-/R(e)), (4.14) 

hJEa 

while the voltage difference is given by 


1/ = ^ (4.15) 
e 

Therefore, the work IV in Eq. fl4.13p is given by 

/•OO 

IV =—— c?er(e)(/L(e) - /R(e)) (4.16) 

e h Jeq 

= c^er(e)(e-/iL)(/L(e)-/R(e)) - der(e)(e -/iR)(/L(e) -/R(e)) (4.17) 

= Jq-Jq, (4.18) 

which results in 

jR 

r,= l-^. (4.19) 

We show in Appendix C that the upper limit of this efficiency is the Carnot efficiency as is 
expected from the theory of the standard heat engine. 

4.1.4 Affinity 

Definition 

Afhnities, or thermodynamic forces are the forces which drive a system in equilibrium out of it, 
such as the difference in the temperature, the chemical potential, the pressure, and so on. We 
dehne the afhnities as follows. We assume that the entropy of the whole system is a function 
of a set of extensive variables {Xk\'- 


^ = ^(Xo,Xi,X2,...). 


(4.20) 
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The entropy production of the whole system is thus given by 




dt 


^ dS dXk 

^ dXk dt 

k 


^ ^ Jk-d-ki 
k 


where we define the affinity and the corresponding ffux by 


Ak = 
Jk = 


dS 

W~k' 

dXk 

dt 


(4.21) 


(4.22) 

(4.23) 


Let us generalize this dehnition of affinities and ffuxes as follows. When the entropy produc¬ 
tion can be expressed as a sum of the products of an extensive flux and an intensive parameter, 


S = ^ JkAk, 

k 


(4.24) 


we call Ak an affinity and Jk the corresponding flux. We show below examples of Ak and Jk- 


Example 1: Affinities corresponding to Jq and 

Let us consider the setup with two terminals shown in Fig. 12.11 As no entropy is generated in 
the wire, the total entropy production is given by the sum of the entropy production of each 
reservoir: 


= /^lQl + /^rQr 
= -I^lJq + I3rJq 

= (4.25) 

which lets us define the affinities A^ = —jdi^ and A^ = of the flows Jq and Jg, respectively. 
Note that the affinities are not a difference of intensive parameters in this case, which might 
be because Jq and Jq are not conserved quantities. 

Example 2: Affinities corresponding to Je and J^ 

We here explain why we chose the affinities of Je and Jn as Eqs. fl3.36p and (13.371) . Using the 
relations Jq = Je — IJ-rJn and Jq = Je — ^^rJn, we can transform the entropy production as 
follows: 


^ = -f^hJq + PrJq 
= —(3l{Je — hrJv) + Pr{Je — /^rJv) 

= {~Pl + (3r)Je + — /^r/^r) Jw 

= A^Je + AcJni (4.26) 

which leads to the definitions of the affinities A^, = —Pr + Pr and Ac = Pr^ir — /^rA^r of the 
flows Je and Jat, respectively. 
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Example 3: Affinities corresponding to Jq or Jq and Jjq 

We then derive the affinities corresponding to Jq and Jn- Using the relation Jq — Jq = 
(Pr — fiL)JN, we can transform the entropy prodnction as follows: 

^ = -/3l + /3r 

= + I3r[Jq ~ (A^R “ ^J'L)JN] 

= {—Pl + I3r)Jq + — f^R)JN 

+ (4.27) 

which leads to Aq = —/3l + /Sr and AJ^ = /Sr(/U.l — /Rr) of the flows Jq and Jat, respectively. In 
the same manner, we can find the affinities corresponding to Jg and Jat: 

S = AlJ^ + A^J^, (4.28) 

where A^ = -/3l + /5r and A^ = /3l(/Rl - /Rr)- 

The notable point here is that the second affinity in either case is proportional to the 
difference in the chemical potential. This was not the case in Example 2. 


Example 4: Affinities corresponding to flows which are superpositions of Je and 

Jn 


We here consider flows which are general snperpositions of Je and Jj^. 
generalized flows Ji and J 2 by 



Let ns define the 
(4.29) 


where 



(4.30) 


with the condition that a, b, c, d are real and det A = ad — be ^ 0. The affinities of these flows 
are given by 

(Ai A2) = {Ah A,)K-\ (4.31) 

becanse we then have 


S - (A N) 

= (A, 

= AiJi + ^2^2- (4.32) 


For Example 1, the matrix K is given by 



as in 



(4.33) 


40 


(4.34) 


while for Example 3, it is 



We thereby confirm 


= - (—/^L + Pr /^L/^L — /5r/^r) ( 

/iL - /iR ^'^-11 

= (-/?L /3r) , 

— (“/5l + /3r PlI^l — PrI^r) 

= {—^L + Pr /5r(/^L — /^r)) • 


(4.35) 

(4.36) 


(4.37) 


(4.38) 


4.1.5 Counting fields 

Let us consider the counting fields corresponding to the generalized flows Ji and J 2 . We 
transform the counting fields Xc and Xh so that the exponent in Eq. (13.281) . 

XcQn + XhQE = {Xh Xc) ’ (4.39) 

may not change. We can achieve it by defining the new counting fields xi, X 2 by 

(xi X2) = {xh Xc)K-\ (4.40) 

Noting that Qi and Q 2 are transformed as Ji and J 2 , we can show that the transformation 
does not change exponent in Eq. (13.281) : 

XcQn + XhQe = {Xh Xc) 

= (X. 

= XiQi + X2Q2- (4.41) 

Let us present the counting fields of several flows; for Jq and Jq we have (I4.34p . and their 
counting fields and are therefore given by 

iX: xf) = ixk Xc)(j i (4.42) 
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(4.43) 


for Jq and Jat, because of (14.361) . the counting fields Xq and Xtv are given by 

(Xq x\j) = {Xh Xc) • 

We can define the counting helds for Jg and Jat in the same manner. 

4.2 The generalized Onsager relations between the heat 
and particle flows 

In this section, we derive the generalized Onsager relations between several pairs of flows. 

4.2.1 General discussion 

Let us dehne the cumulant generating function for Ji and J 2 by 

X 2 , B) = F{xc, Xh, B) = F{hxi + dx 2 , axi + CX 2 ] B), 
with which we can dehne the transport coefficients of Ji and J2 as follows: 

d(iXi)'‘'d(iX2)’“dAl'dAf 
We can prove the following symmetry: 

^JiJ 2 ±(xi, X2; B) = ±FjjJ 2 ±(-xi + iAi, -X 2 + iAx, B), 

where 

BjiJ2±{Xi, X2', B) = Fj^j^{xi-, X2', B) ± Fj^j^{xi-, X2] —B). 

The proof is as follows: 

J 2 ±(xi, X2]B) = F±{hxi + dx2, axi + 0x2] B) 

= F±{xc,Xh]B) 

= ±F±(-Xc + iAc, -Xh + iAh] B) 

= ±F±[6(—Xi + iAi) + d{—X2 + *^ 2 ); a(~Xi + + c(~X2 + *^ 2 ); B] 

= ±^JiJ2±(-Xi + -X 2 + iA 2 -, B), (4.48) 

where we used the symmetry fl3.38p in the third equality. 

From this symmetry fl4.47p . we derive the relations among the transport coefficients using 
the same procedure as in Eqs. fl3.44l) - fl3.46l) : 

aa±(b) = ± e e C') ( 4 . 49 ) 

ni= 0 n 2=0 V V 2/ 

where 

AA±(B} = ± 4‘t’(-S). (4.50) 

There are an inhnite number of choices of Ji and J2, and hence we can derive an inhnite number 
of corresponding Onsager relations. Of course, not all the choices are physically relevant. We 
show below examples of Ji and J2 which have physical meaning. 


(4.44) 

(4.45) 

(4.46) 

(4.47) 
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4.2.2 The generalized Onsager relations between Jg and Jg 

Let us first consider the case Ji = Jq and J 2 = Jq- Because we have 

{Xh Xc) = {Xq xf) K 



= (Xq + Xq -MLXq “ MRXq) , (4-51) 

the cumulant generating function for Jq and Jq is given by 

jR(Xq, xf, B) = F(-/iLXq - hRxf, Xq + xf] B) , (4.52) 

with which we can dehne the transport coefficients of Jq and Jq in the form fl4.45p . The 

symmetry fl4.46p and the relations fl4.49p follow from it. 

4.2.3 The generalized Onsager relations between Jq or Jq and J^ 

We next consider the case Ji = Jq and J 2 = Jn- Because we have 

(Xh Xc) = {Xq Xn) (J 

= (Xq -hLXg + Xn) , (4-53) 

the cumulant generating function for Jq and Jat is given by 

Bj^j^iXQ, Xn, B) = F(-/rLXQ + x\, Xq', B), (4.54) 

with which we can dehne the transport coefficients of Jq and Jn in the form fl4.45p . We 

will specihcally use the notation G for the present case instead of C hereafter for the use 

in Subsection 4.2.4. The symmetry fl4.46p and the relations fl4.49p again follow from fl4.54p . 
We remark that the relations of G reproduce Eqs. fll.l5l) - fll.l7l) by identifying GJq = Gnn, 
Gqi = Gnq, Gio = Gqn, and = Gqq. However, we need to consider the present treatment 
in order to hnd the relations. 

We can do the same for the case Ji = Jq and J 2 = Jn, for which we will use the notation 
M instead of C. We note here that the linear coefficients of M coincide with those of G, but 
they differ in higher orders. 

4.2.4 Application: Nonlinear Seebeck coefficient 

We here show that we can express the nonlinear Seebeck coefficient easily using the transport 
coefficients or We will use Gq’'^^{B) hereafter, but we can use 

in the same way. The advantage to use the coefficients Gq'’^^{B) or MQf^{B) is that the 
corresponding affinities or 4^ and Hjv contain A/i and AT explicitly in contrast to the 
affinities Ac = This enables us to expand thermoelectric coefficients, which are 

usually related to A/r or AT, more easily. 
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The Seebeck effect is a thermoelectric effect in which the voltage difference V = (jUr — jul) /e 
is brought about by the temperature gradient AT = Tr - Tr 0 . The Seebeck coefficient is 
defined by [ 2 ] 


- 



(4.55) 


Let us define the higher-order Seebeck coefficients as follows. Under the condition in which the 
average of the particle flow (Jjv) is zero, the fcth-order Seebeck coefficient is dehned by [371138] 


-V 


SiAT + S2{ATf + ■■■ 

oo 

Y,Sk{AT)\ 

k=l 


(4.56) 


Let us express the second-order Seebeck coefficient with the transport coefficients {B). 
In order to do this, we first express in terms of AAq under the condition {Jn) = 0 , which 


IS 


(^10 „ ^10 
'-^20 /iL 2 ^in .T .T Or, 


{Jn) — G^qA^ + Gq^Aq -h + 


10 /tL /I L I ^-^02 2 _ Q 


2 ! ^ 


Solving this, we obtain two types of the affinity AJ^ as 

■^10 “ ^ii^Q =*= y/ (^10 + y ~ ^2o(2Go?A^ + Gl^^A^ ) 


- 

^N± — 


r<w 

^20 


(4.57) 


(4.58) 


depending on the sign of the square root. Note that we expanded Jn only up to the second 
order of the affinities because the higher-order terms do not affect the result when we calculate 
the second-order Seebeck coefficient. Before expanding fl4.58l) . we have to choose which solution 
we use. We choose A^j^ if G^ > 0 and A^_ if Ggi < ^ following two reasons. One 

reason is that the voltage would not be zero with AT = 0 unless we choose them in this way. A 
finite voltage with no particle flow and no temperature difference is not normal physically. The 
other reason is that the first-order Seebeck coefficient would not coincide with the well-known 
linear Seebeck coefficient [2]. 

We then expand the properly chosen solution with respect to AJq to the second order: 


A^ - 


a 10 /^10/^102 _ o/^lO/^lO/^lO _i_ /^102^10 

01 /(L _ ^ 02^10 ^^ 01 ^ 10^11 + ^01 ^20 2 , 


^10 

^10 


2G 


103 

10 


which turns out to be independent of the choice. Inserting the affinities 

Tl - Tr at 


Aq — —/3l + /5r — 

= /^r(/^L — /Ir) = 


TlTr TlTr’ 

h-L ~ /iR _ ^ 


T 


R 


for Jq and Jn into Eq. fl4.59p . we have 

ey _ GSiGlf - 2G‘;Gi°Gi° + G‘f G'g f AT\\ 

Tr Gig TlTr 201"=’ \TlTr) ' 


(4.59) 


(4.60) 

(4.61) 


(4.62) 
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which is followed by the expansion of the form fl4.56p : 


^10 _ O/^IO/^IO/^IO _i_ ^102^10 

T/ _ '-^01 \rp I '-^02'-^ 10 “I" '-^01 '-^20/A7^12 

^-enG\r ^ 2eT^T^G\f ^ 

From this we obtain the first- and second-order Seebeck coefficients as follows: 

„ ^ 1 GJ» 

' eTLGlg’ 

1 /^lO/^lO^ _ o/^lO/^lO/^lO _i_ 

_ 1 eTQ2Ur]^Q ZUrQ^Urj^gUr^]^ -|- UrQ^ 0^20 

' “ cTl^Tr 


(4.63) 


(4.64) 

(4.65) 


Let us express the second-order Seebeck coefficient fl4.65p only with linear transport coef¬ 
ficients using the relation of the form fl4.49p . For simplicity, let us assume that there is no 
magnetic filed B = 0. We find from fl4.49p that the relations as fl3.47p - fl3.60p is also valid for 

^ki,k2 . 

'^n/2 ■ 


^10 _ ^20 

'-^20 ~ '-^105 


/^ll _ o/^lO /^20 

Orio — Otq^ 

^10 _ ^11 

'-^02 ~ '-^ 01 ) 


^ G}? = 


/^ll I /^20 
^10 '-^ni 


01 


(4.66) 

(4.67) 

(4.68) 


under no magnetic field. Using these relations, we can rewrite the second-order Seebeck coeffi¬ 
cient fl4.65p as follows: 




/^ll/^102 /^lO/^lO/'/^ll I /^20'\ I ^l0^/^20 

^Ol'-^IO ^Ol^lol^lO + ^OlJ + '-^01 ^10 


eTL^TR 


2 G!»“ 


(4.69) 


We thus expressed S 2 only with the linear-transport coefficients. 

In evaluating S 2 using the expression fl4.65p . it may be difficult to measure Gq® and G 20 , 
which are nonlinear coefficients. Using the final expression fl4.69p . however, we can evaluate 
the nonlinear Seebeck coefficient S 2 by measuring the transport coefficients only in the linear- 
response regime. Indeed, it may seem difficult to observe the heat flow. Recently, however, the 
energy flow has been observed at low temperatures in the experiment [35]. We expect that the 
heat flow and its noise will become able to be observed experimentally in the future. 
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Chapter 5 

Summary and Future Works 


In this thesis, we explained the details of the heat flow in mesoscopic one-dimensional transport 
systems and derived the generalized Onsager relations of flows which are superpositions of the 
particle and energy flows. In particular, we derived the relations among the pairs oi Jn and 
Jq, Jn and J^, and and J§. 

We here remark the use of the energy flow and the heat flow. In the research area of 
thermoelectricity, there seems to be a confusion in using the two flows. The concept of the heat 
flow becomes necessary when electrons do work, particularly when we consider a heat engine. 
We thus need to use the heat flow when we consider the work of electrons. 

The relations that we obtained here will be useful in the future. One possibility is to use 
them in calculating the nonlinear coefficients like the nonlinear Seebeck coefficient which we 
obtained in Chapter 4. In order to describe thermoelectric devices which are not in the linear- 
response regime, the conventional thermoelectric coefficients and the hgure of merit, which are 
dehned in the linear-response regime, should be insufficient. We thus need new coefficients 
which describe nonlinear thermoelectric effects. The relations among higher-order cumulants 
will be useful in evaluating the coefficients. 

Another possibility is the efficiency fluctuation. The research of the efficiency fluctuation ap¬ 
peared recently [39114T] . The authors in Ref. [39] proposed that when we consider the efficiency, 
we should replace the conventional second law 


(ARtot) > 0 (5.1) 

with the fluctuation theorem 

(g-AStot) ^ ('5 2) 

where ARtot is the total entropy production of the system. Combining Eq. fl5.2p with Jensen’s 
inequality 

(g-A5tot) > g-{AStot>^ ('5 3) 

we can show that the fluctuation theorem 05.21) includes the second law 05.ip : in other words, the 
fluctuation theorem is a higher entity than the second law. It then necessitates us to consider 
the higher-order cumulants of the entropy production because the quantity contains 

(AS'^qj), . The entropy production is expressed as in Eq. 04.321) by the generalized 

flows Ji and J 2 dehned in Eq. 04.29p . We thus have to consider the higher-order cumulants of 
the generalized hows. Using the generalized Onsager relations which we derived in Chapter 4, 
we may understand the behavior of the efficiency in the nonlinear-response regime. 

The research of nonlinear thermoelectric devices, particularly as a heat engine, is in the 
dawn. We expect that the generalized Onsager relations that we derived in the present thesis 
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will help the development of this research area and manufacturing high-efficiency thermoelectric 
devices using nonlinearity. 
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Appendix A 


Detailed Calculations for Chapter 3 


In this appendix, we show the details of the calculations in Chapter 3. 


A.l Calculation for Eq. (13.291) 


In this section we show the details of calculation to obtain the relation fl3.29p . We start from 
Eq. fl3.28p . which we reproduce here: 


m 

{QNr},{QEr} 


Inserting Eq. fl3.27l) into the above, we have 


Z{{Xcr},{Xhr};B) 




m 

11 '“" cr Q Nr \hrQ Er^ 
r=l 


(A.2) 
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Using Eq. fl3.26p further, we proceed as 


zaXcr},{Xkr};B)= Y 1 

{QiVr}){Qsr} 

m 

X n {■'K"'- - - (Qk“ - Ce”.a)1 

r=l 

z/,A 




r=l 




z/,A 

X (^f |e 


r=l 

m 


3 '^{XcrQ^^\-\-XchQ^r,x) 


^n\^-iHs JJ" ^i{x<^rQ^S^^^^+XchQ'i%) 

r=l 
m 

yinit | ^iHs rr [p-hx cr ^r-\-Xch^r) 


I'^T) 


cr A^r+Xc/i-H'rij |VI/®"^ 


z/,A 


r=l 

m 




r=l 


u,X 

= Tile^^^V^e-^^^V^^po] 

= Tr[e*-^*U2e-*-^"UVoU^] 

= Tr[poU^e*^*U2e-*^*0], (A.3) 


where 


m 

^(Xcr-^r “t“ "KhrHf) 

exp 

r=l 

2 


(A.4) 


We here used the relations W = QjvrV 1^™'*)’ X = Q^r.Al^A^)’ -^r = 

QE^r,ix Hr = Q^r,\ I^^a"); ^^e fact that po and Ul are commutative, and the invari¬ 

ance of the trace under cyclic permutations. 


A.2 Calculation for Eq. (l3.35l) 


In this section, we prove the symmetry of the characteristic function 


Z{{Xcr}, {Xhr}] B) = Z{{-Xcr + iA^r}, {-Xhr + iAhr}', -B). (A.5) 


We first introduce the time-reversal operator 0 [12]. The operator 0 satishes the following 
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properties: 


(A.6) 

(A.7) 

(A.8) 


QiQ = —i, 

{n\0\n') = 

\h) = 0 |n), 

where |n), |n'), and O are bases and an operator in the Hilbert space, respectively. Let us then 
transform the left-hand side of Eq. flA.Sp . We here write the magnetic-held dependence of H 
and V explicitly. Note that H{B) and V{B) are transformed by the time-reversal operator as 
follows: 


eH{B)e = H{-B), (A.9) 

eV{B)e = V{-B). (A.IO) 

We can therefore rewrite the characteristic function as follows: 

Z{{xcr]AXhrhB) = 

= (^1 \n) (A.ll) 

n 

= Y (^1 \h) 

h 

= Y (^1 \n) 

h 

= Y (^1 \n) 

h 

= (A.12) 

where we assumed that the states {|h)} span a complete set if {|n)} do. 

Because Hr and Nr commute with each other for any r, pp dehned in fl3.15l) . or 


PoiB) = 

where z = [ 1^=1 with Ahr = -^r and A^r 

dehned in fl3.30p . or more explicitly 


n 

r=l 


jAhr^r {^)-\-AcrNr 


(A.13) 


—f^rPr, and the operator V 


V {{Xcr}, {Xhr}] ^) = n 


r=l 


-{XcrNr XhrHr{B)) 


(A. 14) 


also commute with each other. We can therefore write down the two terms in flA.12D as follows: 


l^dXcr}, {Xhrh -B)^p'oV{{Xcr}, {XhrY, -B)^ = 


1 


m 


If®''' cr -\-iXcr)Nr-\-{Ah,r~\~iXhT)Hr{ — B) 

r=l 


{Xh.}; -Bf = n 

r=l 


(A.15) 

(A.16) 
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After simple algebras, we can identify them with 


— V{{—Xcr + {“X/ir + —B)^, (A.17) 

z 

{{-Xcr + iAcr], {-Xhr + ^Ar}; -B)^p'^V {{-Xcr + ^^cr}, {-Xhr + iAhr}] -B)\ (A.18) 


respectively. Equation (1A.12P is then rewritten as 


TT[e^^(-B>V{-By poV{-By 

Xcr A j-, Xhr A i-^hr } j ) 


(A.19) 


if pd commutes with H{—B), where in the hrst line we left out the arguments {—Xcr + iAcr} 
and {—Xhr + i^hr} for brevity. 

As is assumed after Eq. fld.lhh . we here used pd = which indeed commutes with H. 

We thereby have the symmetry 


^i{Xcr{j {xhr{] B) — Z(^{ Xcr i-^cr} y { Xhr i-^hr}] B^. (A.20) 
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Appendix B 


The Sommerferd Expansion at Low 
Temperatnres 

B.l The Sommerferd expansion at low temperatures 


In this appendix, we expand the heat flows Jq and Jg using the Sommerfeld expansion at low 
temperatures [13] . 


B.1.1 The Sommerfeld expansion 

We can expand the integral 

poo 

/ deg{e)f{e) 


(B.l) 


in the form 


'Eg 


poo pp. 2rp2 

/ deg{e)f{e) = / deg{e) + -^g\g) + 0{T^), 
’Eg dE g 


6 


(B.2) 


where /(e) = {1 + exp[/3(e — p)]}“^ is the Fermi distribution function and g{e) is an arbitrary 
function which is continuous and infinitely differentiable with respect to e at e = p, while 
diverges no more rapidly than some power of e as e —?■ +cxd. Here, g'{g) denotes the first 
derivative of g{e) with respect to e at e = /i. The expansion (IB.2j) is called the Sommerfeld 
expansion 


B.1.2 The derivation of the Sommerfeld expansion 

Let us derive the Sommerfeld expansion (1B.2I1 . We first obtain the most rough approximation 
of the integral by setting T = 0; 


de g{e)f{e) 


deg{e) 


'Eg 


'Eg 


(B.3) 
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In order to approximate the integral (IB.lj) when the temperature is small but non-zero, we 
subtract the left-hand side from right-hand side of Eq. flB.3j) : 


deg{e)f{e) - / deg{e) 


' Eq 


' Eq 




degie)f{€)- cie^(e)(l -/(e)) 
Jeq 

. _ r d, 




= T 


J Eq 

'dJ±±l£l-T 

e^ + 1 


rfiii^-EQ) gU-Tx) 

ax -- 

. e^ + 1 


(B.4) 

(B.5) 

(B.6) 

(B.7) 


We here used the transformation x = /5(e —/i) in the hrst term and x = — /3(e —/i) in the second 
term of Eq. (IB.61) . Assuming T Tp, where Tp is the Fermi temperature dehned as Tp = ep 
where ep is the Fermi energy, we have 


- Eo) = » 1, (B,8) 

and 

—^1 for X 1, (B.9) 

e* -f 1 

which show that the integrand of the second term in Eq. (IB.71) is exponentially small for large 
values of x. We can thus extend the upper limit of the integral of the second term in Eq. flB.7p 
from /3(e — /i) to inhnity: 


1 e^ + 1 




+ 1 


This lets us expand the integral I in the form 

g{g. + Tx) - g{n-Tx) 


I = T 


dx- 


= 2g\g)T‘^ / dx 


+ 1 
X 

e^TT 


+ 0(T^) 


= %V)t^ + o(t^), 
0 


(B.IO) 


(B.ll) 

(B.12) 

(B.13) 


where we used the Taylor expansion of g, 


g{fi±Tx) 

- gif^) ±9 {f^)Tx + 2 {Txf + ---, 

(B.14) 

and the integral formula. 




f , X 

/ dx -= —. 

Jo e" +1 12 

(B,15) 

Note that 

// N dg{e) 

d. . ■ 

(B.16) 
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B.1.3 Expansion of the conventional heat flow 

We can expand the conventional heat flow using Eq. flB.2l) as follows: 

/•OO 


— 

J n — , 


de{e - /iL)(/L(e) - /R(e)) 


Q h , 

" JEq 

^ roo POO 

de{e - fih)fhie) - t i de{e - /iL)/R(e) 


h 

1 

h 


’Eg 


'Eg 


-ML ^ 2 j ^2 


_J Eg 


de{e - pl) + 


6 


+ 0(Tl^) 


1 

h 


-MR ^2j^^2 


de{e - pl) + 


.J Eg 


6 


+ O(Tr^) 


-^(TR^-TRV^(hK-hO^+o(n 


= r / “ hR)(/L(e) - /R(e)) 


de{e - /nR)/L(e) 


de{e - /iR)/R(e) 


'Eg 


'Eg 


1 

h 


r-ML ^2^12 


de(e - pr) + 


UEg 


6 


+ 0(Tl") 


1 

h 


r-MR 


de{e - pr) + 




6 


+ 0(Tr^) 




(B.17) 

(B,18) 

(B.19) 

(B.20) 

(B.21) 

(B.22) 

(B.23) 

(B.24) 
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Appendix C 

The Upper Bound of The Efficiency of 
The Mesoscopic Heat Engine 


In this appendix, we show that the upper bound of the efficiency dehned by (14.131) is the Carnot 
efficiency and give an example of the transmission coefficient which achieves the upper bound. 


C.l The upper bound of the efficiency of the mesoscopic 
heat engine 

C.1.1 General upper bound of the efficiency 

Let and denote the entropy productions in the left and right reservoirs, respectively. 
Using the equalities dQ = TdS = dU — /idiV, we can relate these entropy productions to the 
conventional heat flows as 


J^ = Ul- /iriVL = 

Jg = f/R - /TrAr = Tr^r. 


These relations let us transform the efficiency fl4.13p to the form 


V 


IV_ 

tL _ jR 
/L 

-Tl^l - Tr^r 
Tl^l + Tr^r 


(C.l) 

(C.2) 


(C.3) 

(C.4) 

(C.5) 

(C.6) 
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We here impose the condition that the net entropy production S' = S'l + (Sr is positive, that is, 
S'r > — S'l, which gives 


V = 

< 


Tl^l + Tr^i 


R 


- Tr^l 


= 1 - 7^ = he, 
-l-L 


(C.7) 

(C.8) 

(C.9) 


where r]c is the Carnot efficiency. We can achieve the equality if and only if S'l = S'r, that is, 

^ = 0 . 


C.1.2 Example of the transmission coefficient which gives the Carnot 
efficiency 

The Landauer-Biittiker formalism with the transmission coefficient gives 

7 = f / der(e)(/L(e) -/R(e)), (C.IO) 

^ JEg 

1 

Jq = T c?e^(e)(e-/^L)(/L(e)-/R(e)). (C.ll) 

" JEg 

We thus express the efficiency fIC.dp in the form 

IV 

V=-fr (C.12) 

_ (/iR - /in) /g der(e)(/L(e) - /r( 6)) 
c?er(e)(e -/iL)(/L(e) -/R(e)) ■ 

We now know that we can achieve the Carnot efficiency when the total entropy production 
of the system is zero. We thus express the total entropy production with the Landauer-Biittiker 
formalism: 


^ + ^R 

_ Qi^ Or 

“ Tl Tr 

tL tR 

= __Q _L iR. 
Tl Tr 


deT{e) 


I Eg 

nOO 




/^R 


(/L(e) - /R(e)) 


Tl Tr 

der(e){[log/L(e) - log(l - /L(e))] - [log/R(e) - log(l - /R(e))]} 


J Eg 

X (hie) - fnie)) 


c?er(e)(/L(e) - /R(e)) log 


'Eg 


/L(e)(l - /R(e)) 

./R(e)(l - /L(e)) 


(C.14) 

(C.15) 

(C.16) 

(C.17) 

(C.18) 

(C.19) 












Let us assume r(e) > 0, which is physically reasonable. We can show that the integrand in 
Eq. fIC.lQj) is always non-negative as follows. When /L(e) > fnie), 


which leads to 


log 


/L(e)(l-/R(e)) 


_/R(e)(l - /L(e)) 

When < /R(e), on the other hand, we can show 

7L(e)(l -/R(e)) 


log 


_/R(e)(l -/L(e))_ 


fhie)fRie) 

(C.20) 

/L(e)/R(e) 

(C.21) 

- fhie)), 

(C.22) 

> 0. 

(C.23) 

> 0 

(C.24) 


similarly. Using these inequalities and r(e) > 0, we can show that the integrand is always 
non-negative: 

7L(e)(l -/R(e))' 


^(e)(/L(e) - /R(e))log 


> 0. 


(C.25) 


_/R(e)(l-/L(e))_ 

Let us then consider the condition for S = 0. We easily hnd the following condition; for 
each value of e, r(e) = 0 or fhi^) — fnie) = 0. If r(e) = 0 for any e or if fhi^) — fnie) = 0 
for any e, however, the transport would not happen, that is, the condition is trivial. We thus 
have to hnd a nontrivial condition in which r(e) ^ 0 and fhi^) — fni^) = 0 at an energy; if we 
demanded fhi^) = /R(e) at two energies, they would be equal at any energy. 

In order to do this, let us set the transmission function m 


r(e) = 6{e - ej. 


(C.26) 


where Cc is a constant. This transmission function is not a sufficient condition but a necessary 
condition for us to obtain S' = 0. Let us also set the condition 


fhiec) = fniec), 


which gives the value of Cc as follows: 


fhiec) = fniec) 

1 




1 _l_ e^Lhc-ML) 1 -p g/^ahc-MR) 

f^hiec — hr) = f^Riec — /xr) 

^LhR — ^RhL 


er. = 


(C.27) 

(C.28) 

(C.29) 

(C.30) 

(C.31) 


Tl-Tr 

Substituting Eq. flC.26p into Eq. flC.13p . we obtain the following expression of the efficiency: 

(hR - hr) JZ c?er(e)(/L(e) - fnie)) 


r] = 


IZg deTie)ie - /iL)(/L(e) - fnie)) 
(hL - hR)(/L(ec) - /R(ec)) 


hL)(/L(ec) - fniec)) 


/^L “ /^R 
/^L 


(C.32) 

(C.33) 

(C.34) 
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Substituting Eq. flC.31j) into Eq. flC.34j) . we indeed achieve the Carnot efficiency: 


/^L — /^R 

r] = - 

— /^L 

/^L — 

_ ,, 

7 l-Tr 

_ ^L/^R ~ ^L/^L ~ ^R/^R + ^R/^L 
^L/^R — ^L/iL 

_ ^l(/^r — /^l) — T^if^R — hr) 
^l(/^r — /^l) 



Note that when rj = rjf.^ the electric current and the heat current vanish: 

e 1'°° 

^ = T c?er(e)(/L(e) - /R(e)) 

" JEg 

= ^(/l(6.)-/r( 6,)) (•.•r(6)=<5(6-6,)) 

= 0 (■.•/L(ec) =/R(ec)), 

1 

'^Q = T c?er(e)(e-/iL)(/L(e)-/R(e)) 

= ^(ec - /^L)(/L(ec) - /R(ec)) C.' r(e) = 6{e - ej) 

= 0 (•.• /L(ec) = /R(ec)), 

which gives the vanishing power IV. This is the same as the standard heat engine; 
cycle produces zero power. 


(C.35) 

(C.36) 

(C.37) 

(C.38) 

(C.39) 

(C.40) 

(C.41) 

(C.42) 

(C.43) 

(C.44) 

(C.45) 

the Carnot 
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